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Abstract. Using the auxiliary field method, we give an analytical expression for 
the eigenenergies of a system composed of two non-relativistic particles interacting 
via a potential of type \J d^r'^ + h. This situation is usual in the case of hybrid 
mesons in which the quark-antiquark pair evolves in an excited gluonic field. 
Asymptotic expressions are proposed and the approximate results are compared 
to the exact ones. It is shown that the accuracy is excellent. 
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1. Introduction 



The auxiliary field method (AFM) has been recently developed and used to compute 
analytical approximate relations for bound state eigenvalues of the Schrodinger 
equation. Formulas very accurate are obtained for power-law potentials [Ij, sums 
of two power-law potentials [2] and exponential potentials [Sj. More recently it has 
been shown that, although obtained in very different ways, the AFM method and the 
envelope theory [3] are completely equivalent [5]- Let us note that this equivalence 
leads to a deeper understanding of both frameworks. 

The aim of this report is to give an analytical expression of the eigenenergies of 
the Schrodinger equation with the potential 

V{r) = Va\^ + b, (1) 

using the AFM. Such an interaction has a strong interest in hadronic physics, in 
particular for hybrid mesons in which the quark-antiquark pair evolves in an excited 
gluonic field. The applications of the results obtained here for potential ([T]) to such 
systems will be given elsewhere [6]. 



2. Eigenenergies 

2.1. Analytical expression 

Let us follow the general procedure of the AFM |lj . Our goal is to find approximate 
expressions for the eigenvalues of the Hamiltonian 

H^^ + Vir), (2) 
2m 

where m is the reduced mass of the particles and V{r) is given by ([1]). We first choose 
an auxiliary function P(r) = r^; the auxiliary field v is then defined by 

V'(r) 

y ^ Kir) ^ = ^^^=. (3) 

For the moment v is an operator, and ([3]) can be inverted to give r as a function of v. 
r = I{v). Explicitly 



The AFM needs the definition of a Hamiltonian H{v) = p'^/{2m) + vPir) + 
V{I{v)) — vP{I{v)). In our particular case, 

H{.) = f + ur-+^^ + %. (5) 

If we choose the auxiliary field in order to extremize H: 5H / 5v\^^o — 0, then the 
value of this Hamiltonian for such an extremum is precisely the original Hamiltonian: 
II{v) = H. Instead of considering the auxiliary field as an operator, let us consider it 
as a real number. In this case, the eigenenergies of H are exactly known for all (n, I) 
quantum numbers: 



m w a 

where, as usual, N — 2n + I + 3/2 is the principal quantum number of the state 



= + - + (6) 
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The philosophy of the AFM is very similar to a mean field procedure. We first 
seek the value of the auxiliary field which minimizes the energy, dE/dv\^^^g, and 
consider that the value E{vq) is a good approximation of the exact eigenvalue. It is 
useful to use the new variable 

and to define the parameter 

The minimization condition is concerned now with the a:o quantity and results from 
the fourth order reduced equation 

44 - 8x0 - 3y = 0. (9) 

The solution of this equation can be obtained by standard algebraic techniques. It 
looks like 

xo = G{Y) = l^v(Y) + yA{V{Y))-y^-V{Y), (10) 

with 

v{Y) = (2 + + ya) ^'^ - y (2 + ^4 + r^) ^'^ . (11) 

Substituting this value into the expression of E{vq) leads to the analytical form of the 
searched eigenenergies, namely 

E{vq) = 

nr r _ 11 

(12) 




G{Y)_ 

The problem is entirely solved. 

As it is shown in [5], the same formula would be obtained for the choice 
P(r) = sgn(A)r'^ with A > —2, but with different forms for the quantity N . With the 
choice A = 2 made above, N = 2n + l^ 3/2. In this case, using results from ,5 , it can 
be shown that formula (fT2|) gives an upper bound of the exact result. For A = — 1, 
N — n + l + 1 and the formula gives a lower bound. The qualities of these bounds are 
examined below. The form of N is not exactly analytically known for other values of 
A. 

An approximate simpler form of (|12p avoiding the complicated G function and 
giving the lowest order exact results in both limits y ^ and y ^ 00 for finite value 
of 6, is given by 




^('^O) ~\h;T7\ \/3y + (3 X 22/3 _ + ^ 1 (13) 



where rj is an arbitrary parameter. A very good approximation is obtained for 77 
around 1: for a fixed value of 6, the relative error between and is below 2%. 
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2.2. Asymptotic expansions 

At long range, the potential ([1]) behaves as the linear potential ar. This asymptotic 
behavior is equivalent to the limit 6-^0. In this case F — > but h/{?>Y) 
{2a^N'^/my^^/16. Moreover G{Y) 2^/^ when Y ^ 0. Reporting these conditions 
in the value i?(fo) given by (fT2|) . one obtains the asymptotic behavior 

Ei^or^K-Y^N^/'- (14) 

This is precisely what is expected for a pure linear potential ar (see i.e. [1]). Very 
accurate estimation of the exact eigenvalues can then be obtained by using, for 
instance, iV = (7r/\/3)n + I + \/37r/4 [HE]. 

Another interesting asymptotic expression is the limit — )■ oo, or equivalently 
a — > 0. It is easy to check that, in this limit, the potential ([T]) reduces to 

V{r)^Vb+^^. (15) 

Thus the potential is just equivalent to a harmonic oscillator plus a constant term. 
Consequently, the exact energies are given by 

E* = " (2n + l + l)+Vb. (16) 



When a 0, then Y oo and it is easy to see that G{Y) « (3^/4)^/'*. For this 
hmit, ini) reduces to IH]) for the choice N = 2n + l + 3/2. In this case also, the AFM 
leads to the good asymptotic behavior. 

3. Scaling laws 

Let us denote by E*{m, a, b) the exact eigenvalue of the Hamiltonian The scaling 
properties of the Schrodinger equation (see [2]) allow to express it in term of the 
eigenenergy of a reduced equation. More precisely 

where e* (/?) is the exact eigenvalue of the dimensionless Hamiltonian 

HiP) = ^ + v/^^. (18) 

Switching to the AFM approximation, one can verify that the approximate energy 
E{m, a, b) as given by p2)) satisfies the same scaling law ([T7|) as the exact energy, the 
reduced approximate energy e(/3) being still given by (|12p in which E = e, b — f3 and 

Y = ^ (19) 

In consequence, to test the quality of the approximation it is sufficient to make 
comparisons between e(/3) and £*(/?)• This is the subject of the next section. 

The Fourier transform of the Hamiltonian H = + V a^r^ + b is a spinless 
Salpeter Hamiltonian with a harmonic potential 

Hq^ = UJ^/I^1~^P + ar^, (20) 
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where the following substitutions have been made [7] 



4a , , Vb 



mauj 



M=—, a=^—. (21) 



In order to get a relevant equation for particles of mass M, uj must be set to 1 or 2. So 
the results obtained here can also be used to study the spectra of the Hamiltonian ([^0]) . 
Such a task will be developed in another work where the AFM will be applied to 
relativistic Hamiltonians [8]. 



4. Comparison to exact results 

As remarked previously, the AFM cannot give strong constraints on the dependence 
of A'' in terms of {n,l). In particular, had we chosen P{r) — sgn(A)r'^, the better 
choice for A^ would have been A^ = A{X)n + I + C(A), with the quantities A[\) and 
C(A) given in [1]. The square root potential ensures a smooth transition from a linear 
form (A = 1 but in this case we have only approximate expressions) to a quadratic 
form (A = 2 and in this case the values are exact) as (3 increases from to oo. 

It is thus natural to suppose a smooth dependence of these coefficients on the only 
relevant parameter of the problem, namely [3. Therefore we choose, for A^ appearing 
in Y through (|19p . an expression in the form 

N = A{l3)n + l + C{l3). (22) 

From the results of [1], it is expected that A{Q) « 7^/^/3 « 1.814, C(0) « \/37r/4 « 
1.360, lim;3^oo A{I3) = 2 and Yivap-^oo C{l3) = 3/2. 

The procedure we adopt is based on the following points: 

• We calculate the exact values e*;(/3) for < n < rimax, < Z < Zmax and for a 
given set of P values. This program is achieved using a very powerful method 
known as the Lagrange mesh method (described in detail in [9]). For our purpose, 
we consider that rimax — knax = 4 is a good choice. For any calculated value, we 
have an accuracy better than 10~^. 

• We calculate the approximate values e„i(/3) using (fT2|l with Y given by (fT9| in 
which A^ is deduced from (|22p. for the same set of (3 values. Building the x-square 

-, "max 'max 

we request the coefficients A{(3) and C(/3) of A^ to minimize this function. The 
obtained values are represented by black dots in figures [TJ 

• In order to obtain functions which are as simple as possible, continuous in (3, and 
which reproduce at best the above calculated values, we choose hyperbolic forms 
and require a best fit on the set of the sample. Explicitly, we find 

_ 8/3+102 _ 30/3 + 53 

4/3 + 57' ^^^^" 20/3 + 39- ^^^> 
These integers are rounded numbers whose magnitude is chosen in order to not 
exceed too much 100. The corresponding values are plotted as continuous curves 
in figures [TJ They have been constrained to exhibit the right behavior A —^ 2 and 
C 3/2 for very large values of /3. Formulas ^ give A{0) = 102/57 « 1.789 
and C(0) = 53/39 « 1.359. Both values are very close to the theoretical numbers 
given above. 
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Figure 1. Best values of the coeflScients A{(3) and C{0) to parameterize tlie 
eigenvalues of Hamiltonian l|18p : numerical fit with II23II (dots); functions l|24p 
(solid line). 



Since our results are exact for /3 — s- oo, one has obviously x = in this limit. The 
error is maximal for small values of /? but, over the whole range of (3 values, the results 
given by our analytical expression can be considered as excellent. Just to exhibit a 
quantitative comparison, we report in table [1] and in figure [2] the exact e* ; (/?) and 
approximate e„;(/3) values obtained for /3 = 1, a value for which the corresponding 
potential is neither linear nor harmonic. As can be seen, our approximate expressions 
are better than 1% for any value of n and / quantum numbers. Such a good description 
is general and valid whatever the parameter [3 chosen. 



Table 1. Comparison between the exact values e*;(/3) (2nd line) and analytical 
approximate expressions (/3) for the eigenvalue of Hamiltonian l)18|l with /3 = 1. 
For each set (n, I), the exact result is obtained by numerical integration. 3rd line: 
approximate results are given by ||12|I with l|19| l. \22) and l|24|l : 1st line: upper 
bounds obtained with N = 2n + I + 3/2; 4th line: lower bounds obtained with 
N = n + l + l. 





I = 


/ = 1 


I = 2) 


/ = 3 


I = 4 


n = 


1.94926 


2.49495 


2.99541 


3.46197 


3.90193 




1.91247 


2.45074 


2.94841 


3.41419 


3.85430 




1.89549 


2.44621 


2.95032 


3.41969 


3.86189 




1.65395 


2.22870 


2.75000 


3.23240 


3.68492 


n = 1 


2.99541 


3.46197 


3.90193 


4.32027 


4.72059 




2.89556 


3.34652 


3.77899 


4.19405 


4.59335 




2.85420 


3.32970 


3.77678 


4.20097 


4.60620 




2.22870 


2.75000 


3.23240 


3.68492 


4.11355 


n = 2 


3.90193 


4.32027 


4.72059 


5.10556 


5.47723 




3.74112 


4.14232 


4.53310 


4.91307 


5.28251 




3.69078 


4.11913 


4.52783 


4.91998 


5.29790 




2.75000 


3.23240 


3.68492 


4.11355 


4.52250 


n = 3 


4.72059 


5.10556 


5.47723 


5.83725 


6.18692 




4.50374 


4.87138 


5.23246 


5.58628 


5.93264 




4.44883 


4.84403 


5.22459 


5.59242 


5.94903 




3.23240 


3.68492 


4.11355 


4.52250 


4.91485 


n = 4 


5.47723 


5.83725 


6.18692 


6.52732 


6.85935 




5.20859 


5.55148 


5.88996 


6.22329 


6.55111 




5.15078 


5.52098 


5.87970 


6.22821 


6.56756 




3.68492 


4.11355 


4.52250 


4.91485 


5.29295 
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Figure 2. Spectra e„i of Hamiltonian 1181 1 with /3 = 1. For each value of n, 
the eigenvalues are presented for I varying from to 4. Diamonds: exact results 
obtained by numerical integration. Circles: approximate results given by I I12I I 
with I I19I I. II22I I and II24I I. The error bars extend from lower bounds obtained with 
N = n + l + lto upper bounds obtained with N = 2n + I + 3/2. 



The upper bounds obtained with P{r) — are far better than the lower bounds 
computed with P{r) = — 1/r. This is expected since the potential \/ a^r^ + b is closer 
to a harmonic interaction than to a Coulomb one. Better lower bounds could be 
obtained with P{r) = r. But, the exact form of N is not known for this potential, 
except for I = for which N can be expressed in term of zeros of the Airy function. 
With the approximate form N — (7r/\/3)n + I + -\/37r/4 [HE], we have checked that 
results obtained are good but the variational character cannot be guaranteed. 

5. Conclusions 

In this paper, we propose an analytical approximate expression for the eigenenergies of 
a Schrodinger equation for two non-rclativistic particles interacting via a potential of 
type V a^r^ + b. This situation corresponds to the case of a hybrid meson in which the 
quark-antiquark pair evolves in an excited gluonic field [HI . We give the corresponding 
expressions for any value of the parameters a and b and for any values of the radial n 
and orbital / quantum numbers. Thanks to a Fourier transform, the energy spectrum 
wc find can also describe a relativistic one-body or two-body Hamiltonian with a 
harmonic potential. 

The scaling laws properties are shown to be fulfilled exactly by these approximate 
expressions; moreover the limiting cases 6 — > (linear potential) and 6 — > oo (harmonic 
potential) reduce to the exact solutions. 

The approximate analytical results are compared to the exact ones. It is shown 
that for any values of the parameters and for a whole range of quantum numbers the 
obtained accuracy is excellent. The formulas we get are expected to play an important 
role in the identification of hybrid mesons ^ . 
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